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Abstract 

The next-to-leading order cross section for the inclusive production of charm quarks 
in 77 collisions is calculated as a function of the transverse momentum pj- and the 
rapidity y in approaches using massive or massless charm quarks. For the direct 
cross section we derive the massless limit from the massive theory with the result 
that this limit differs from the massless version with MS factorization by finite 
corrections. Subtracting or adding these corrections allows us to compare the two 
approaches on equal footing. We establish massless and massive versions with 3 
and 4 initial flavours which are shown to approach the massless approximations 
very fast with increasing px- With these results we calculate the inclusive D*^ 
cross section in 77 collisions using realistic evolved fragmentation functions with 
appropriate factorization scales and compare with recent data for da/dpx from 
three LEP collaborations after single- and double-resolved contributions have been 
added. 
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1 Introduction 



Recently the three LEP collaborations, ALEPH [|I|], L3 and OPAL [|] have presented 
data for inclusive D*^ production in two-photon collisions at e^e^ center-of-mass energies 
close to \/S = 189 GeV. Besides the total cross section a for 7 + 7 ^ + X also the 
differential cross sections with respect to the D* transverse momentum, da/dpjn, and the 
pseudo-rapidity, da/dr], have been measured. 

In 77 collisions, where both photons are on-shell, each of the two photons can behave 
as either a point-like or a hadron-like (or resolved) particle. Therefore one distinguishes 
in such collisions three production channels called direct (both particles interact point- 
like), single-resolved (one 7 is point-like, the other hadron-like) and double-resolved (both 
photons are hadron-like). The resolved contributions require the knowledge of parton 
densities in the photon, whereas the production via two direct photons is free of such 
non-perturbative input. The different channels mix in higher orders of perturbation the- 
ory and thus the distinction between the direct and resolved contributions becomes scale 
and scheme dependent: only the sum of all three contributions is a physical cross section. 

The transverse momentum distribution da/dpx for the inclusive production of D*"^ mesons 
is characterized by two distinct scales, the mass m of the charmed quark and the transverse 
momentum of the D* or charm-quark. Depending on the ratio px/^i two different 
approaches for next-to-leading order (NLO) calculations in perturbative QCD have been 
used for a comparison with the experimental data 0, 0, In the so-called massless 
scheme, the charm quark is considered as an active flavour in the photon 0. In these 
calculations, it is assumed that four flavours q = u, d, s and c are present in the pho- 
ton, described by corresponding distribution functions, and all quarks are taken to be 
massless. The charm-quark is also an in-going parton originating from the photon in 
the case of the resolved contributions to the cross section and it fragments into the D* 
meson similarly as the produced u, d, s quarks and the gluon g. The predictions of this 
approach are expected to be valid in the region of large transverse momenta pr ^ rn. 
In this scheme, calculations for the small px region are not reliable. The cross section 
diverges in the limit px ^ and the total cross section can not be predicted. Following 
the usage in deep inelastic charm production we shall refer to this as the zero-mass 
(ZM) 4-flo,vour scheme (it corresponds to the "massless charm scheme" introduced in Ref. 
in connection with charm production in 7p collisions). 

The other scheme, in which cross sections for 77 — >■ D* X have been calculated [Q], is the 
so-called massive charm scheme 0. In the massive charm scheme the number of active 
massless flavours in the initial state for the resolved contribution is equal to n/ = 3 and 
the charm quark is assumed to be massive. The massive c quark appears only in the final 
state. In this scheme the charm mass, m ^ ^qcd, acts as a cutoff for the initial and final 
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state collinear singularities, and sets the scale for the perturbative calculations. The cross 
section factorizes into a partonic hard scattering cross section multiplied by light quark 
and gluon densities in the case of the resolved contributions. In leading order (LO), the 
direct production is described by the partonic reaction 7 + 7 — > c + c while the resolved 
contributions involve the channels 7 + (7 — > c + c (single resolved) and g + g — > c + c and 
g + g ^ c + c (double resolved), where q are the light (massless) quarks q = u, d, and s. 
This approach has the advantage that not only the various distributions, like in rapidity 
and/or transverse momentum, can be predicted in the full range of pt but also the total 
cross section. 

One might expect that the massive approach is reasonable only in those kinematical re- 
gions where the mass m and any other characteristic scale like pt are approximately of 
the same magnitude and significantly larger than Aqcd- Under these circumstances the 
charm mass can be used to set the renormalization scale entering the quark-gluon coupling 
as as well as the factorization scale needed to evaluate the quark and gluon densities of 
the photon. In NLO, terms oc as ln(p|./m^) arise from collinear emission of the gluon 
by charmed quarks at large transverse momentum or from almost collinear branching of 
photons or gluons into cc pairs. These terms are not expected to affect the total pro- 
duction rates, but they might spoil the convergence of the perturbation series and cause 
large scale dependencies of the NLO result at pt S> m.Q In the massive approach the 
prediction of differential cross sections is thus limited to a rather small range of pt — ttl. 
Nevertheless, predictions of this approach have been compared to data up to pt — 10 
GeVi. 

The proper procedure for p-r ^ m is to absorb potentially large logarithms into distribu- 
tion and fragmentation functions where they can be resummed by virtue of the Altarelli- 
Parisi equations. To implement this procedure one needs a charm contribution in the 
photon parton distributions (PDF) and a fragmentation function (FF) for the transition 
c ^ D* (or any other charmed meson or baryon). Logarithms oc ln(M^/m^) defined with 
the factorization scale M are absorbed into these distribution and fragmentation functions 
and remaining terms oc ln(p|./M^) are of order 0(1) for the appropriate choice M ~ pt- 
For sufficiently large pr the cross sections calculated with this finite charm mass method 
including four active fiavours must approach the results described earlier with the zero- 
mass 4-flavour scheme 0. Since m 7^ 0, it is reliable also for intermediate pt > ttl, where 
most of the experimental data have been obtained so far, rather than only for pr ^ m. 
We shall call this scheme with a massive charm quark, but with terms proportional to 
as ln(M^/m^) subtracted, the NLO 4-fl(ivour scheme. The scheme described above, where 
this subtraction is not performed and only three light active fiavours are considered will 
be called the NLO 3- flavour scheme. This latter heavy quark mass approach is also of- 
ten referred to as the fixed-fiavour number scheme in connection with heavy quark 

^ Similar potentially large terms cx as lii{Q^ /m^), Q being the photon virtuality, appear in the calcu- 
lation of charm electro-production cross sections at next-to-leading order Q . 
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electro-production. In the following we shall not use the term "massive charm scheme" 
any more, since it is ambiguous. Instead we shall refer either to the NLO 3-flavour or the 
NLO 4-flavour approach. 



The relation of these three approaches for the process 77 D*X has not been investi- 
gated in detail yet. In particular, the relation between the ZM 4-flavour and the NLO 
4-flavour scheme, e.g. the question, in which range of pt values the ZM 4-flavour scheme 
is a good approximation to the NLO 4-flavour scheme, has not been studied so far. Also, 
it is of interest to know, for which px the NLO 3-flavour and the NLO 4-flavour scheme 
produce approximately the same results. So far only differential cross sections calculated 
in the ZM 4-flavour and in the NLO 3-flavour scheme have been compared 0. It was 
found that the two approaches differ in the definitions and relative contributions of the 
direct and resolved terms, but essentially agree in their sum. However, in this comparison 
0], the ZM 4-flavour result was modified in so far as terms containing ln(M^/m^) have 
been taken into account not only at order 0{as), but resumming them with the Altarelli- 
Parisi equations in the so-called perturbative fragmentation function approach, in which 
initial conditions for the FF's were taken from perturbation theory This made it 

difficult to pin down the finite charm-mass effects present in the NLO 3-flavour scheme. 



It is the purpose of this work to fill this gap and present a comparison of results obtained 
in the ZM and the m 7^ 4-flavour schemes. We will identify terms in the massive theory 
surviving in the limit m — > which are not present in the ZM approach where the quark 
mass is put to zero from the beginning. These terms describe final-state interactions and 
can be interpreted as a perturbative fragmentation function describing the transition from 
massless to massive charm quarks. Only after correcting for this difference by subtracting 
the final-state interaction terms from the massive theory (or adding them to the massless 
theory) one can expect that both theories approach each other in the large-p^- limit. This 
way, the m = and the m 7^ theories can be considered on the same footing and a com- 
parison will be sensible. In addition we shall compare also to the NLO 3-flavour scheme in 
order to find out the difference at small and intermediate pt- We shall concentrate in this 
comparison on the direct cross section for 77 — Z)* X, since this is the dominant part. A 
similar comparison for the single-resolved and the double-resolved cross section is left to a 
later study. The merging of the massive fixed-order approach and the resummed massless 
fragmentation approach has been investigated in Ref. for hadron-hadron and photon- 



hadron scattering. In principle one can infer from these results the corrections due to the 
finite charm mass for the single- and double-resolved cross sections. This work, however, 
as well as the calculations for 77 reactions in Ref. 0, is based on the perturbative FF's 
pi| approach, which is too restrictive for our purpose. 

The outline of our work is as follows. In section 2, we shortly describe the formulae which 
are used to calculate the cross section for 7 + 7— i^c/c + X with m 7^ . We derive from 
these cross sections the limit m ^ and compare it with the results of the zero-mass 
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theory which can be found in the hterature. This defines the necessary subtractions in the 
massive theory so that it approaches the massless theory in the hmit of large pt ■ Section 
3 contains the numerical results for the comparison of the two theoretical approaches 
based on different choices for the scales at which finite initial and final state terms are 
subtracted or absorbed into non-perturbative PDF's of the photon or FF's of the charm 
quark. In this section we also present comparisons to the 3-flavour scheme. After adding 
single- and double-resolved contributions we compare the results for the m ^ and the 
m = 4-flavour schemes to recent experimental data from LEP II. Our conclusions are 
summarized in section 5. 



2 Decomposition of the LO and NLO Differential 
Cross Section 

2.1 Leading-Order Cross Section 

We first consider the process 

7(Pi) + 7(P2) ^ c{ps) + c{p4) + [g{k)] (1) 

where pi, i = 1, 2, 3, 4 and k denote the momenta of the two incoming photons and the 
outgoing c, c quarks and a possible gluon (in square brackets), which is present when we 
consider the NLO corrections. Below we will describe the procedure needed to obtain 
differential cross sections for D* production in 77 scattering. We have the following 
invariants 

s = {pi+ P2)^, t = T -w? = {pi- ps^ - rn^, u = U - = {p2 - Pzf - (2) 
and 

S2 = 'S'2 - = (pi +P2- Pzf -ni^^s + t + u. (3) 
It is customary to define the dimensionless variables 

V ^ 1 + -, w ^ (4) 

s s+t ^ ' 

so that 

t = — v), u — —svw^ S2 = sv{\ — vo). (5) 
The leading-order cross section is 

— -!p- = c(s)5(l-w) - + - + 4 4 (6) 



where 

27iNca^ei 

c s = . (7) 

s 

Nc is the number of quark colours and Cc is the electric charge of the charm quark, 
Cc = 2/3. From (j^) the finite charm mass corrections are clearly visible. In the next 
section we shall show numerical results explicitly. 



2.2 The Next-to-Leading-Order Cross Section 

The NLO corrections consist of two parts, the virtual corrections to 7 + 7— s^c + c and 
the gluonic bremsstrahlung contributions 7 + 7 — ^ c + c + g. These NLO corrections have 

ni PI . Only in Ref . explicit formulae for the 



been calculated by several groups [|T^ 
separate contributions due to one-loop diagrams and due to bremsstrahlung contributions 
are given in a form which allows us to derive the massless limit (m ^ 0). The results in 



Ref. []T3| are subdivided into three parts, the vertex plus self-energy cross section daysE, 
the virtual box cross section dasox and the gluon bremsstrahlung cross section dasr, so 
that the NLO single-inclusive differential cross section da/dvdw is decomposed as follows 



da 



NLO 



daysE _|_ da-Box _^ da-sr 



dvdw dvdw dvdw dvdw 



The three parts in (H) have according to Ref. |jT3[ the following structure 



da 



VSE 



dvdw 



C{s) 
4 

- In 



S{l-w)\2Ai 4 



8-6 



T 

(2 - Lia ( ^ 
t 

' T 



m 



+A2 In f ^ ) + ^3 



Li2 



T 



C2 



+ + (t 



u 



(9) 



and 



daBox CpO-s dai^o 2rn? — s 
dvdw 7T dvdw s(3 



{-2 In X In /3 + 2Li2(-a;) - 2Li2(a;) - 8(2} 



+ ^^5 1 - w){ -85i ^ Inxln — 

4 y sp Vm^ 

Q / \ f ^ 

+2— ( Inx 41n(l + x) - Inx - 41n ( — 



-t 



+ 4Li2(-x) + 2C2 



T 



+2B3ln' X + 4-^\nx + AB^ln { — ] +8^6 In — +4^7(2 + 453 
p \m^J \m^J 

+ {t^u)\. (10) 
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In (|T0|) 



and in 



(3 = y^l — Aw?/ s , X 



1 + 



and (0) C2 = and the normalisation 



C{s) 



C[S] 



27r 



(12) 



The quantities Ai and Bi are functions of m^, s, t and m. They are given in appendix B 



of Ref. |jT3| and will not be repeated here. The contributions (P) and ([T0|) contain also in- 
frared divergent terms proportional to (2e = 4 — n) in dimensional regularization with 
dimension n. They are omitted since they cancel against terms in the bremsstrahlung 
cross section. How these terms are distributed in rfuvsE and (icTBox can be inferred from 
Ref. [1131. 



The last term in (gj) looks more complicated. According to Ref. |13] it can be written in 
the following form: 



dvdw 



i SVS2 ( S2{S + U\ 

'145^'^^^^ + 



250 



AS. 



S2{s + U 



1 ^2 ^ 

In —63 + 



4^2 



rn?[s + uY 



64 



+1565 + JsCs + /gCg + /loCio + /laCis + /isCis + AeCiG + (t ^ m)) 



1 



ei 



(l-«;)+452 
+52/12612 + ^2/14614 + it 



2S2, T 
— ^ In ■ 

y 



U-y 4S2 

rCe H ^67 + 52/11611 



r + f/ + y 



u 



-2 



27r dvdw s/3 
-4/? 



(2m — s) 41nxln — - + 21nx 



Lio 



+ In a; 



+ 2sP 



2 In 



(13) 



Here we used y = \/ (t + uf — Am?s and the coefficients ci, 62, 63, 64, 65, ce, 67, es, eg, 
eio, eii, ei2, 613, 614, eis and eig are again functions of the invariants s, t, m, S2 and of rn?. 
I5 to /16 are integrals over angles, which have been evaluated in Ref. |]T3| and are written 
down in appendix C of this reference. 



The cross sections daysE and dasox are proportional to 5(1 — w). The bremsstrahlung 
cross section da^i contains terms proportional to 5(1 — w) and to Other con- 

tributions are finite for w — > 1, as long as m 7^ 0. In the limit m — they give rise 
to additional terms proportional to 5(1 — w), (13;;;) > as well as terms proportional to 
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A second approach where the mass of the charm quark is neglected from the beginning 



was worked out in Ref. |T5[ and later confirmed in Ref. ||T6[. With m = 0, coUinear 
singularities appear. They are regularized by dimensional regularization. In order to 
understand the differences of the two approaches, we have evaluated (|), ([T0|) and ([l3|) 
in the limit m —>■ 0. Special care must be exercised in order to recover all the terms 
proportional to 6{1 — w), and ■ We write the result in a form which has 



been introduced in the calculation for massless quarks in Ref. |T^. This will allow us to 
identify the terms which come in addition to the massless theory in the MS factorization 
scheme |jT6|. The LO cross section for m = has the simple form 



lim ^'^^'^ = c(s)(5(l — w)ro(t') with tq^v) = 1 . (14) 

dvdw 1 — V V 

The decomposition of the NLO cross section in the limit m ^ has the form 

hm — — ; — = Cl + Cl m — \d[l — w) 

■m^o dvdw \ s J 

( ~ , / 1 \ /ln(l -w)\ 

+ C2 + C2 In — + C3 ' ' ^ 



s J\l — wJ+ \ 1 — w 
+C5 In t> + Cg ln(l — vw) + C7 ln(l — v + vw) + Cg ln(l — v) 

+C9 In u; + Cio ln(l - w) \ Cn + \ cx\ + Cn j In ■ 



2 
m 



s 

l-v 



\n{l-v + vw) \nw Hjz^) 

+C12 : hci3- hci4— . (15) 

I — w I — w 1 — w 



The coefficients ci, ... C14, ci, C2, cn, and cn are either functions of v alone or of v and 
w. We obtained the following results for them: 



where 



Cl 



C{s) (^ln\l - v) 
+ ln^f 
+ f4C2 



2 1 

+ 



3 2 
+ --3 

1 — V V 

7 



V 1 — V 
+ hiv 



\n{l-v) 
3 



l-v 



1 



+ 7^ + 2 



2(1 -v) 2v 



To{v) + Aci 



(16) 



Aci = C{s) ( -2 In^ w - 2 In w + 2) to{v 



(17) 



~Ci = C{s) i^-2\nv-^yo{v); 
8 



(18) 



where 



where 



C2 = C(s) (21n^;-^)ro(^;) + Ac2, (19) 



Ac2 = C(s)(-41n^;-2)ro(^;); (20) 



~C2 = -2C{s)to{v) ■ (21) 



C3 = 2C(s)to(v) + Ac3 , (22) 



Ac3 = -AC{s)tq{v) ; (23) 



c _ C7(s)p^*^^ + ^2w(l-t;) 

y(l — vu')^ [1 — vwY [1 — v){l — vw) 1 — vw 

2v 2 , 3 2 2v ^ , 

Av-\ \ \-2w -\- Avw 



1 — V V 

where 



1 — V + vw V vw w w 

+ — + Ac5 , (24) 



^, , / 2v 2f2ty 4 2f 4 4f \ ,^ , 

Acs = C s + + + ; 25 

\1 — f I — V w w vw I — V + vw I 



C6 = 2C{s) [ ,^ V (2v^ -2-3v^ + 5v) + (v^w - 2v^) ] ; 

\v[l — V)W ^ ^ 1 — V ^ ' ) 



(26) 



C7 = -C(.) i + — ) ; (27) 

\1 — 1) + vw vw J 



(28) 



cg = '^-!^^''^\ {v'^w'^ + v'^w - (1 - v)) ; 

v[l — v)w ^ ' 



= C5 + ^^^^ , f4u^w - Iv^vp- -2v^ -2 
w{l — V) ^ 



^, J 2v 8v(l - v) 7v 

cii = C{s)' 



'^l — V + vw {1 — vwy {1 — vwy 

T 1 

-71 ^1 (2-7y + Sy2^^+4v-2 + - + — 

[l—v)[l — vw)^ ' V 1 — 

/ o 8 2 \ 2v 1 1 \ ^ 

-w{v-2 + - + + + - — + Acn , 

V V 1 — V J w w zvw I 



Acn = C s + - + - w(l + v) 

\1 — V + vw 1 — V 1 — V 

2 V 2 \ 

+ 1 ; 

w w vw J 



^, , / 2v{\-v) 2v V 2vil 



{1 — vwY [1 — vwY {l — v){l — vw) 1 

V ^ 2 1 V v'^w 2w\ 
+ 4v-\ 2vw + ; 



1 — V V vw W 1 — V V J 



^, , / 2v V 1 + (1 — v)^ v'^w 
cii = C{s) ■ + ^ + 



1 — V -\- vw 1 — V vw 1 — V I ^ 



1 + v^ 

ci3 = 2C{s) 



v{l — v) ' 
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ci4 = 2C{s) 



1+11 
v{l - 



(37) 



By comparing with appendix A of Ref. |TB| we obtain agreement if we set in (|TBp the 
factorization scales M| = M|. = vn? and put Aci = Ac2 = Acs = Acs = Acn = in 
the formulae above. This means, that the m — >■ limit of the NLO correction as derived 
in Ref. [|13| agrees with Ref. |]I6| with initial factorization scale M| = w? and final state 



factorization scale M% 



m 



and all Acj terms put to zero. Denoting the ZM NLO result 



of Ref. |jT6[ with dazu/dydw, we can write 



with 



da 



da 



NLO 



da- 



ZM 



dvdw dvdw 



{Mi = Mf = m) + 



da- 



FSI 



dvdw 



(3^ 



FBI 



dvdw 



Aci(5(l -w) + Ac2 




Acglnt; + Acii . (39) 



In Ref. 16] terms containing initial and final state singularities have been subtracted using 



different factorization scales M/ and M^?, respectively. Therefore, a separation of terms 
proportional to ln(M|/s) (coefficient Cn) and ln(M|./s) (coefficient cn) was possible. 
Starting from a calculation with m 7^ 0, this separation can, of course, not be obtained - 

]) we have split 



both terms appear with the logarithm ln(m /s). In our formulas (^ 



this sum in two terms so that they agree with cn and cn in Ref. [|T^. In Ref. |jT6[ the MS 



factorization scheme is defined with the spin averaging for the incoming photons taken 
to be 1/2(1 — e) and not equal to 1/2 as it is sometimes done. From our calculation we 
conclude that the zero-mass limit of the massive theory automatically yields the result for 
the MS scheme with n-dimensional photon spin averaging. By comparing with Ref. [|r^ 
we note that Ci and C2 originate from the final-state logarithmic terms. The initial-state 
logarithms have no factors proportional to 5(1 — w) or 1/(1 — u')+ and appear only in Cn, 
i.e. in the non-singular contribution for w ^ 1. There are also no terms proportional to 
ln(/i^/s), where /i is the renormalization scale, since as is not renormalized in first order. 

c/c + X in the limit m ^ 

only 



Ml 



m 



Our derivation shows that the massive theory for 7 -|- 7 
approaches the massless theory in MS factorization with scales Mf 
if the terms Aci, Ac2, Acn are removed. This result is to be expected since the 
regularization of coUinear singularities with a mass parameter m does not necessarily 
give identical results as with dimensional regularization and m = from the beginning. 
The difference is a sum of finite terms, which have first been derived by calculating the 
cross section for e"*" -|- e~ c/c + X in the two theories [|10|- In Ref. it was shown 
that the additional finite terms can be interpreted as a partonic fragmentation function 
d'^{x, iio) for the transition from massless to massive charm quarks, c(m = 0) — > c(m 7^ 0). 
The massless theory with MS factorization has to be folded with d'^{x, /io) in order to 
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obtain the cross section in the zero-mass hmit of the massive theory. This partonic or 
perturbative FF in 0{as) is 



27r 



1 — X 



In 



pi 



- 21n(l - x) - 2 



(40) 



One can show that the extra terms da-psi/ dvdw can be recovered in the framework of the 
massless theory in Ref. |T^ by convoluting d'^{x,fio = m) as a FF of the charm quark 
with the massless LO cross section^: 



dvdw 



hm 



dahoipT = Pt/x) 
dvdw 



dl (x, fio = rn) 



(41) 



This means that the Acj terms represent final-state interactions which survive in the 
massive theory even in the limit m — > 0. Therefore we have labeled their sum with the 
subscript "FSI". They are non-singular, i.e. not proportional to ln(m^/s). 



In Ref. the expression (|40| ) was used as a starting condition at the scale = m for 
the calculation of perturbative FF's at an arbitrary scale via the usual Altarelli-Parisi 
evolution equations in the massless approach. Our derivation shows that we obtain the 
same finite final-state interaction terms as in e^e~ c/c X. One expects that they 
appear in any process for which cross sections in the limit of the massive theory and the 
massless theory with common factorization are compared in NLO. We remark that such 
finite terms do not appear in connection with initial state factorization, since the distri- 
bution of quarks in the photon has identical non-singular terms for the two cases where 
one considers the massless limit of the massive theory and the theory with massless quarks. 



The final result for the massless limit of the massive theory, as given by (^) with the 
coefficients written down in (|^) to (^), can be used in two ways. On the one hand, 
we can compare the cross section of the original massive theory with the massless theory 
by calculating the cross section from the full formula (|D together with (H), (p!OD and 
(P^) with the limit result in (|T3p. In this way we can establish how the additional terms 
proportional to m? and possibly m? Inm^, i.e. all terms which vanish for m? 0, depend 
on the kinematic variables. In particular, we can find out, in which kinematical region, 
the massless limit in (p!5| ) is a good approximation to the massive theory. Compared to the 
massless cross section as derived in Ref. [|l6l the cross section in (p!5[) includes the finite 



final-state interaction terms Aci, Ac2, Acn, and is computed for the factorization 
scales Mf = Mp = m^. On the other hand, it is clear that the massless approximation 
(p!5|) is relevant only for ^ m?. In this case is the large scale, so that m? is not a 
good choice for the factorization scale. The appropriate choice is M| = M|, = [rn? + p^) 
with ^ = 0(1). Thus it makes much more sense to compare the two theories, massive and 



^The explicit form for this convolution and prescriptions for the rescaling of the kinematic variables v 
and w can be found in llSll. 
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massless, for this choice of the factorization scales. Then the large logarithms oc ln(s/m^) 
are removed and absorbed in the PDF of the photon or in the FF of the D* meson. These 
PDF's and FF's are in NLO usually constructed in the MS factorization scheme which 
we shall assume in the following. In order to establish the effect of terms proportional 
to m? in this more realistic case, we subtract the finite final-state interaction terms Aci, 
Ac2, Acii from the massive theory and change the factorization scale in the massive 
theory to M| = Mp = {w? + p^) by additionally subtracting 



2 

m 



Ci In -772"5(1 — w) 



f da \ 

2 1 2 2 

~ TTl 1 TTl ~ TTl 



from the massive cross section above or equally from the massless cross section in (|T^ 
which then is identical to the result in Ref. [|16|. Of course, the difference of the mas 



sive and the massless approximation will be the same as in the first case with m as 
factorization scale, only the relative amount will change. It is clear that the two cross 
sections, massive or massless, for M| = Mp = ^ (m^ + pp) must be supplemented with 
a realistic choice for the c —>■ D* FF and with the additional single- and double-resolved 
cross sections having a charm contribution in the photon PDF. We remark that the latter 
approach makes sense only in the NLO 4-flavour scheme whereas the first approach with 
Mf = Mp = is meaningful only in the NLO 3-flavour scheme. 



3 Numerical Results 



We start with presenting results for the Pt distribution of the cross section for 

e+(p+) + e-(p_) -.c(p3)+X. (43) 

For this, we have to fold the cross section given in the previous section with the quasi-real 
photon spectrum given in the Weizsacker- Williams approximation by 



ZTT X 



+2(1 -x) 



X) 



In 



E^e^Jl - x)2 + mW 



mlx'^ 



mix 



E'^ei{l - xf + mlx"^ 



X 



(44) 



Here E is the e+ (e^) beam energy, 9^ the maximal angle under which the outgoing 
electrons (positrons) are tagged and x = E^/E the fraction of the beam energy entering 
the 77 cross section. Fixing the transverse momentum pt of the charm quark and its 



13 



rapidity y, one can perform the integration over the energy fractions of the two photons 
Xi = Ei/E, i = 1, 2, with the result: 

da 2pT dv /•! dw ,. / VW\ . f l-V\ da 



dpTdy S Jvw f (1 — v) Jvw/v w \ vw J \ 1 — v J dvdw 
where 



1 

=xiX2S 

(45) 

y=^-f-^e-y, W = ^f--^ey, S = AE\ (46) 

In the calculation with massive charm, the charm mass is kept non-zero everywhere, also 
in the definition of kinematic variables and phase space limits. In particular, y is the 
rapidity defined by u = \ For the comparison between various versions of the 

NLO cross section and their massless limits we consider the pt distribution da/dydpx as a 
function of pt for y = 0. Since the cross section as a function of y for fixed pt is maximal 
at ?/ = we expect a similar behavior for the cross section da/dpx after integration over y. 

In all the following comparisons the charm mass has the value m = 1.5 GeV. The total 
energy is = 2E = 189 GeV and the maximal angle in ( ^i)) is 9^ = 0.033 in accordance 
with the choice in the OPAL experiment [0. The normalization of the cross section, as 
given by (^5]) and (|^) (or (|14D) in the LO case is changed by multiplication with the 
branching ratio BR{c — > D*) = 0.235. This value for BR is the average of the five exper- 
imental measurements collected and discussed in Ref. [0 {BR = 0.235 ± 0.007(±0. 007)). 
A more recent value for this branching ratio based on measurements of the four LEP and 



the SLD experiments is BR = 0.241 ±0.008 |T8|. The calculated cross sections are for the 



sum of c and c production. Including the branching ratio for c D* in the cross section 
will make it easier to compare the results for the cross sections without explicit FF for 
c — >■ D* with the cross sections with FF included at the end of this section. We should 
stress that the figures of this subsection do not include the resolved contributions so that 
a comparison with experimental data is not yet reasonable. 



We start with the LO cross section calculated from (^Sf ) with da /dvdw given in (P) and 
its limit for m = in ([T^) . The result for the two cross sections da/dydpx for m 7^ and 
m = 0, respectively, is shown in Fig. |I]a. For large pt the two cross sections approach each 
other. Their ratio, da{m)/dydpT/da{m = 0)/dydpT, i.e. the massive over the massless 
cross section, as a function of pr, is plotted in Fig. [l|b. As is seen in this figure, the 
deviation of the massless from the massive cross section is appreciable at small pt- The 
massless cross section increases strongly for pt ^ 0, whereas the massive cross section 
increases only moderately with decreasing pt, making the ratio smaller than one. For 
Pt > 5 GeV the massive cross section differs by less than 10 % from the massless cross 
section. In Fig. |I|b the ratio crosses the line 0.9 at px = 4.9 GeV. 

In the next two figures. Fig. ^ and Fig. ^ we show the same comparison including the 
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Figure 1: (a) pr distribution da/dydpx for e^e~ 77 c/c + X —> D*^ + X (direct 
contribution, BR{c D*^) = 0.235) at LO for y = in the massless (full line) and the 
massive calculation (dashed line) with m = 1.5 GeV; (b) ratio of the massive and massless 
calculations. 



full NLO corrections. The value for as is calculated from the two-loop formula with 
A(Ar^=4) = 328 MeV which corresponds to as{mz) = 0.1181 in accordance with Ref. ||T9| . 

The renormalization scale is always /ir = \Jpt + in all the results of this and the 
following section. 

In Fig. ^ the cross section da/dydpx at ?/ = is shown in the NLO 3-flavour scheme. 
For m 7^ this is the full cross section with da/dvdw in ( ^5)) calculated from (|8|-|13D. In 
the limit m — > we have added the final-state-interaction coefficients Aq and changed 
the factorization scales to Mj = Mp = m. Only with these changes, we can expect that 
massive and massless cross sections approach each other for large px- This is indeed the 
case as seen in Fig. ^,b. If we compare with the LO cross section in Fig. |l]a we notice 
that the NLO cross section is smaller than the LO one, in particular for increasing pt, 
i.e. Pt > 3 GeV. This reduction of the NLO 3-flavour cross section is due to the choice 
of the factorization scale Mj = Mp = m, which, of course, is the appropriate scale in the 
case of fixed three fiavours. This generates negative terms in the NLO corrections which 
increase in absolute value with increasing pp. These large correction terms proportional 
to ln(m^/s) (see (|T^)) make the perturbative cross section more and more unreliable, the 
larger pp is, and call for a subtraction of these terms by choosing a factorization scale 
determined by the large scale pt- The subtraction terms originate predominantly from 
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Figure 2: (a) pt distribution da /dydpx in the NLO 3-flavour scheme (FSI-coefficients Aq 
added in the massless calculation, direct contribution including BR{c D*^) = 0.235) 
for y = in the massless (full line) and the massive calculation (dashed line) for Mj = 
Mp = m; (b) ratio of the massive and massless calculations. 



the collinear singularities in the final state and, only to a smaller extent, from the coUinear 
singularities of the initial state (see (^) for the separation of the two parts). We remark 
that we used the same definition of in both schemes, i.e. we took Nf = 4 and the same 
value for A(7v^.=4) in the formula for as, although, for consistency with the choice of three 
flavours, we should have taken Nf = 3 and the corresponding value for A(Ar^=3). This 
prescription will make it easier to compare with the results in the NLO 4-flavour scheme. 
The correct choice with Nf = 3 in would have decreased the cross section shown in Fig. 
^ by a few per cent. The ratio of the massive to the massless cross section in the NLO 
3-flavour scheme is shown in Fig. as a function of Pt- The approach towards unity 
with increasing px is slower than for the LO ratio in Fig. |I]b. This is presumably caused 
by additional terms oc m?\n.m? which are not present in LO. Now the ratio is > 0.9 for 
Pt > 7.3 GeV, i.e. in the NLO 3-fiavour scheme, larger values of pt are required if the 
massless cross section should be a reasonable approximation. Towards smaller px the 
massless cross section overestimates the full massive cross section quite strongly as was 
already noticed in Ref . . Ai px = 2 GeV this overestimation amounts to approximately 
85%. 

The cross sections in the NLO 4-flavour scheme are plotted in Fig. |^a, again for m 7^ and 
m = 0. For m 7^ the cross section differs from the one shown in Fig. Qa by subtraction 
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Figure 3: (a) pr distribution da/dydpr in the NLO 4-flavour scheme (FSI-coefficients 
Acj = in the massless calculation not resummed and subtracted in the massive calcula- 
tion, direct contribution including BR{c — D*^) = 0.235) for y = in the massless (full 

line) and the massive calculation (dashed line) for Mj 
the massive and massless calculations. 



Pj< + m?; (b) ratio of 



of the FSI-coefficients and by changing the factorization scales to M/ = Mp = \Jpt 

i.e. we subtracted from the massive cross section (p|-p!3D the terms Aq and the logarithms 



The massless cross section is obtained from (O, 16-37) with Aq = 0, i.e. it is iden- 



tical to the results of Ref. []TB|. We observe that the NLO 4- flavour cross section is larger 
than the NLO 3-fiavour cross section in Fig. |]a. About half of this difference is due to the 
enlarged factorization scale; the other half is due to the subtraction of the FSI-coefficients 
Acj. The increase is most effective at the larger pt values. Without additional fragmen- 
tation effects due to the c — > D* transition the NLO 4-flavour scheme would thus lead to 
larger cross sections than the NLO 3-flavour scheme. The ratio of the massive and the 
massless cross section in the NLO 4-fiavour scheme is presented in Fig. §3. It is clear that 
the difference between these two cross sections is equal to the one in the NLO 3-fiavour 
scheme shown in Fig. |]a, but the ratio changes slightly. Now, the ratio is larger than 0.9 
at Pt = 6.5 GeV, i.e. at somewhat smaller px than in the NLO 3-flavour scheme (Fig. |^). 

It is clear that the results presented up to this point can not be compared to experimental 
data as long as no FF for the transition c D* is included. This deficiency is amended 
now for the cross section daiy = 0)/dydpT shown in Fig. ^a and Fig. ^ with two different 
choices of the FF. The cross sections in Fig. |a, again for the massive and the massless 
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version, are the cross sections in the NLO 4-flavour scheme augmented with the so called 
perturbative fragmentation function (pFF) The pFF is evolved via the Altarelli- 

Parisi equations to the appropriate scale with the initial condition at scale /xq equal to 

D:{x,fio)=S{l-x)+d:{x,fio) (47) 

where (i^(x,/io) is given in (|40|) with the choice fiQ = m as initial scale. This prescription 
amounts to taking into account FSI terms not only at order 0{as), but to resum them 
by virtue of the Altarelli-Parisi equations. The normalization of the cross sections again 
includes the experimental branching ratio as in the previous figures. Such a FF can be 
considered as a perturbative model for the real FF which we can use instead of a cal- 
culation based on non-perturbative methods. Alternatively, a FF can be extracted from 
experimental data for the production of D* mesons obtained from other measurements. 
The latter approach will be discussed below. In Fig. ^ the perturbative FF is applied 
equally to the massive and massless cross sections. Through the Altarelli-Parisi evolution 
the terms oc as{fi)lia{MF/ fio) are resummed, which leads to an appreciable reduction of 
the cross sections at large pt- Comparing to the cross sections in Fig. ^ the reduction 
amounts to a factor of more than two at px = 12 GeV. In the small-p^ region the influ- 
ence of the fragmentation is largely reduced. The cross section with the pFF based on the 
boundary condition (|47| ) was used in the massless approximation for a comparison with 
the massive NLO 3-flavour theory in Ref. Comparing with Fig. ^ we see that the two 
cross sections are approximately equal at pt — 5 GeV and the NLO 3-fiavour scheme has 
a much larger cross section than the ZM 4-fiavour scheme at pt = 12 GeV in agreement 
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with the resuhs in Fig. ^ of Ref. 0. The difference at larger pt is due to the choice of 
scheme and the effect of the evolution of the FF in the ZM 4-flavour scheme. Of course, 
as in the previous figures the ZM cross section approaches the massive one for large pt- 
The ratio of the two cross sections with m ^ and m = in the NLO 4-flavour (pFF) 
approach is exhibited in Fig. ^ showing that the ratio goes to one for large pt- The ratio 
is larger than 0.9 for p^ > 5.6 GeV, i.e. the crossing point occurs at a smaller px than in 
the NLO 4-flavour cross section without FF (see Fig. |]b). 



The approach based on perturbative fragmentation functions is, however, not sufficient 
to describe the fragmentation c — > D*. For example, in Ref. [^] it was shown that in 
order to account for the inclusive production of D* mesons in e~^e~ annihilation at various 
center-of-mass energies an additional non-perturbative component is needed. This means 
that the fragmentation of charm quarks into D* mesons cannot be calculated in perturba- 
tion theory. A non-perturbative component, which is not known theoretically and has to 
be determined from other data, is always needed. Hence, it is more appropriate, to give 
up the perturbative component of the FF input altogether and to describe the c —>■ D* 
transition entirely by a non-perturbative FF, as it is done for the fragmentation of u, d 
and s quarks into light mesons. This approach has been followed in Ref. ||^. There a 
non-perturbative FF for c D* was constructed by fitting data on e+ -|- ^ D* + X 



at the Z mass. For the fits two data sets from ALEPH |^ and OPAL ||2^ have been 
used. Since these two data sets are not compatible with each other, two separate fits 
have been performed. In the following we shall employ the results obtained from the fit 
to the OPAL data 0. LEP data are best suited for such fits since they have a better 
accuracy as compared to lower energy data (see Ref. for details). However, with LEP 
data there is the additional complication that one has to separate two different sources 
of charmed hadrons: in hadronic Z decays charmed hadrons are expected to be produced 
either directly through the hadronization of charmed quarks in the process Z —> cc or 
via weak decays of B hadrons produced in Z ^ bb with approximately equal rate. With 



the data from Refs. |21, 22] it was possible to disentangle these two components and to 



construct the FF for c — >■ D* by determining the FF at the starting scale /io = 2m. At 
this initial value the charm-quark FF was assumed to have the Peterson form, which is 
particularly suitable to describe EE's that peak at large x. The Peterson FF depends only 
on two parameters, the normalization factor and the shape parameter e. According 
to Ref. H the fit to the OPAL data gave in NLO the values = 0.267 and e = 0.116 
resulting in a branching ratio BR{c D*) = 0.238 at = mz in good agreement with 
the average value from Refs. [1^, IT! 



Based on the purely non-perturbative FF of Ref. |Q (OPAL set at NLO) we have cal- 
culated da/drjdpT at r/ = for massive and massless quarks as a function of pr- tj is 
the pseudo-rapidity which was used also in the analysis of the experimental data ^ |^ 
and we identify the pseudo-rapidity of the D* with the rapidity of the charm quark. The 
resulting cross sections, denoted NL0-4(BKK), are shown for 2 GeV < < 12 GeV in 



19 



1 1 1 r 



d^tPb/GeV] 
NL0-4(BKK) 




J I I L 



1.2 



1 - 



: 0.8 - 



: 0.6 - 



i 0.4 - 



0.2 



10 12 

PT[GeV] 



1 


1 1 1 1 




da{m) jda{m=Q) 




dydpx ' dydpx 




y = 


1 


NL0-4(BKK) 
1 1 1 1 



6 8 10 12 



Figure 5: As in Fig. ^ but including c D* fragmentation according to Ref. [§] (set 
OPAL NLO). 



Fig. ^a. The cross section is smaller than the one in Fig. This results essentially from 
the fact that the non-perturbative FF peaks at a value x < 1, whereas the perturbative 
FF is dominated by the contribution proportional to 5(1 — x). For a given transverse mo- 
mentum pt of the observed D*, smaller x in the FF probes larger px = Pr/x of the charm 
quark in the underlying hard scattering process where the cross section is smaller. The 
approach of the massive to the massless cross section is even faster than in the previous 
studies. This is seen in Fig. where the ratio is plotted. The ratio exceeds 0.9 already 
at pt = 4:.7 GeV. This stronger reduction of mass corrections is also due to the different 
x-dependence of the fragmentation functions since mass terms proportional to oc iv? /p^ 
decrease faster with increasing pt than the massless cross section. 

In Fig. 1^ we show how the NLO results in Figs. 0, ^ ^ and ^ are related to each other. 
For this purpose we introduce a common normalization by dividing the results of all four, 
massless and massive, cross sections by the massless LO cross section (full curve in Fig. 
|l]a). The full lines in Fig. ^represent the massless calculations for NLO-4, NLO-3, NLO- 
4(pFF) and NL0-4(BKK). The dashed curves correspond to the massive cross sections. 
Compared to the LO massless cross section, the NLO corrections in the NLO 4-flavour 
scheme produce an increase of the cross section between 40 % (low pt, massless) and 
20 % (large pt)- The increase at low pt is compensated by mass effects and at large pr 
by including the pFF so that the net effect of NLO contributions in the range 3 GeV 
< Pt < 12 GeV is a reduction between 30 % and 50 %. Most of the reduction of the 
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Figure 6: pt distributions normalized to the result of the LO massless calculation (full 
curve in Fig. |l|a), da/dai^oim = 0) at ?/ = 0. Full lines: massless calculations, dashed 
lines: massive calculation. From top to bottom: NLO 4-flavour scheme with Mj = Mp = 
+ m? (NLO-4); NLO 3-flavour scheme with Mi = Mp = m (NLO-3); NLO 4-flavour 

scheme with M/ = Mp = ^p"^ + including perturbative fragmentation (NL0-4(pFF)); 
and NLO 4-flavour scheme with non-perturbative fragmentation (NL0-4(BKK)). 

NLO 3-flavour result compared to the NLO-4 curve originates from the different choice 
of the factorization scale. At low pp, around ~ 2 GeV, the massive cross sections are 
approximately equal to each other. The inclusion of a fragmentation function leads to 
a further reduction of the cross section which is particularly strong for the case of the 
non-perturbative FF. 

From this section we conclude that the corrections due the finite charm mass are ap- 
preciable at low pp and are small, i.e. below 10 % for pp > 5 GeV. Furthermore, the 
fragmentation corrections due to the evolution of the perturbative FF are also large, but 
now in the large-p^ region. For the non-perturbative FF there is also a reduction at smaller 
Pp due to the fact that in this case the hard cross section is probed at considerably larger 
transverse momenta of the charm quark. The approach of the massive theory towards the 
massless approximation is always very fast. It is strongest for the theory with a realistic, 
i.e. non-perturbative FF. For this case, dX pp = 2 GeV, the neglected mass terms lead to 
an overestimation of the massive result of not more than 40 %. It is also apparent, that the 
pure NLO 3-fiavour approach, i.e. without any additional non-perturbative fragmentation 
corrections can give a trustworthy prediction only for rather small pt, i-e. for pp <?> GeV. 
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At higher px the effect of evolution is essentiaL But even at small pt-, we find a reduction 
by a factor 2 due to the non-perturbative FF, i.e. even at small px the NLO 4-flavour 
approach with perturbative FF does not lead to a realistic prediction for the cross section. 



4 Comparison with LEPII Data 

Experimental data for the differential cross section da/drjdpT integrated over some fixed 
f] region come from the three LEP collaborations ALEPH, L3 and OPAL. The ALEPH 
data |l|] represent da/dpT integrated over —1.5 < r/ < 1.5 and are averaged over the 
LEPII runs with in the interval 183 GeV < < 189 GeV. The most recent L3 
data are integrated over |?7| < 1.4 and averaged over 183 GeV < < 209 GeV. 
The more recent OPAL data (second Ref. in [Q) are integrated over |?7| < 1.5 with the 
luminosity averaged = 193 GeV (183 GeV < < 202 GeV). We compare the three 
data sets of Refs. |2^ with our calculation at \fS = 189 GeV. Besides the differing 
a/S* regions there are the slightly different regions of the r] integration and possibly differ- 
ent anti-tagging conditions for the outgoing electrons and positrons. We disregard these 
differences and put the data points of all three experiments in one plot. We have checked 
that the different l^^l max chosen by L3 on one side and ALEPH and OPAL on the other 
side change the cross section by at most 7 % which is much below the measurement errors 
of the data points. Also the different values of y/S are not expected to change the cross 
section by more than the corresponding experimental errors. As an example we checked 
that increasing from 189 to 193 GeV increases da/dpr by 2 (3)% at pr = 2 (12) 
GeV. The influence of different anti-tagging conditions could not be investigated because 
of insufficient information in the corresponding references. 



For the comparison with the experimental data we need also the contributions of the 
single-resolved and double-resolved channels. They are available only in the massless ap- 
proximation. Therefore we shall add the direct cross section as calculated in the previous 
section (in the massive and massless version) and the single-resolved and double-resolved 
cross sections in the massless approximation. For this we use the codes developed in Ref. 
for the calculation of the single-resolved and double-resolved cross section for the pro- 
cess 77 — > 71^ X. They were also used in Ref. |2^ and in the first reference of |^ for the 



process 'jp D X. In all cases we shall use the BKK fragmentation function as de- 



scribed in the previous section. We choose as factorization scales M/ = Mp = ^ yPr + 
with ^ = 2. This allows us to calculate da/dpx down to smaller pt than previously since 
the starting scale of the non-perturbative FF of Ref. was put equal to 2m. 

The partition of da/drjdpT into direct and resolved contributions, integrated over t] in 
the region —1.5 < rj < 1.5 is shown in Fig. ^ Here, DD, DR and RR denote the direct, 
single-resolved and double-resolved cross sections in the massless approximation and in 
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Figure 7: Direct and resolved contributions to the pt distribution da/dpT after integration 

over \r]\ < 1.5 at NLO with Mj = Mp = 2\Jp\, + m? including BKK-fragmentation. Up- 
per long-dashed hne: direct, dotted line: single-resolved, lower short-dashed line: double- 
resolved contributions, full line: sum of all contributions. 

NLO, respectively. At p-r = 2 GeV, these three contributions amount to approximately 
38 %, 21 % and 41 % of the total sum. At pr = 12 GeV the relative contributions are 
71%, 14% and 15%, respectively. This means that the resolved cross sections DR and 
RR decrease with increasing pt much faster than the direct component DD. The separate 
contributions depend strongly on the factorization scales: for M/ = Mp = \Jpp + m?, i.e. 
with ^ = 1, the direct contribution is larger by roughly 10 % for all pp-, the single-resolved 
contribution becomes steeper (with an increase of 38 % at = 2 GeV and 3 % at = 12 
GeV), whereas the double-resolved contribution is much reduced (by 60 % at pt = 2 GeV 
and 30% at pp = 12 GeV). We give these numbers for completeness only, although the 
scale dependence of the separate contributions is unphysical. The scale dependence of the 
physically observable cross section is weak: for the sum of all contributions these changes 
combine to a 13 % decrease at pp = 2 GeV and a 1 % increase at pp = 12 GeV. The 
single-resolved cross section has an appreciable component originating from the charm 
PDF in the photon. It is equal to 15 % at = 2 GeV and equal to 61 % at p^ = 12 GeV. 
The double-resolved cross section is due entirely to the charm component in the photon. 
In Fig. and the following figure the photon PDF is taken from Ref. [^, which has an 
explicit charm contribution. 

The sum of all three contributions, DD, DR and RR, is compared to experimental data 
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Figure 8: px distribution da/dpT after integration over |?7| < 1.5 in the NLO 4-flavour 
scheme with Mj = Mp = 2^Jpj< + m? including BKK- fragmentation compared to LEP 
data [H, 0, Full line: massless calculation, dashed hne: massive calculation. Single- 
and double-resolved contributions are included using photon PDFs of Ref. I2H . 



in Fig. ||. The full curve is the cross section in the massless approximation as in Fig. ^ In 
the dashed curve the direct massless cross section is replaced by the direct cross section 
with massive quarks, i.e. NL0-4(BKK) of the previous section, except for the change of 
factorization scales. The resolved components are as in Fig. 0. The experimental data 
points shown at pt values between 1.5 GeV and 8.5 GeV are from ALEPH |l[], L3 
and OPAL 0]. The overall agreement between the theoretical prediction (dashed curve) 
and the experimental data is quite good although the data points in the medium p^ range 
lie slightly above the theoretical curve. Even if a finite charm mass correction for the 
DR and RR contributions would be included, for which we expect a reduction of the 
theoretical prediction by approximately 15 % at pt = 2 GeV and less at higher p^, the 
overall agreement for pt > 2 GeV would hardly change. The data from OPAL 0] and L3 
23| have been compared already with the predictions of the massless theory, which was 
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obtained on the basis of the work in Ref . P . Compared to these results we neglected the 
fragmentation of the gluon g ^ D*. We checked that these contributions are very small 
at small pr and amount to a contribution of approximately —1 % at the largest pt, which 
is negligible compared to the errors of the experimental data. This is due to the fact that 
the gluon enters only at NLO and that the g —>■ D* FF is small compared to the dominant 
charm FF. A similarly small effect results from a change of the charm mass. Choosing 
m = 1.3 GeV increases the cross section by 1 % in the considered pt range between 2 
GeV and 12 GeV. 



5 Summary and Conclusions 

In this work we compared two approaches, the massive and massless schemes, for the 
calculation of inclusive charm production with 3 and 4 initial flavours. As a first step, the 
cross section da/dpT is calculated for the direct component of the 77 reaction in NLO. We 
found that the massless limit of the massive cross section differs from the massless theory 
with MS factorization by finite terms which are non-singular for m — > 0. A first sensi- 
ble comparison of the results of the massless theory with the massive 3-fiavour approach 
can be performed after adding these finite terms to the massless theory. The additional 
terms can be interpreted as the order 0{as) coefficient of the perturbative fragmenta- 
tion function describing final-state interactions of the transition from massless to massive 
charm quarks. It is clear that a resummation of these terms by using a perturbative frag- 
mentation function is not sufficient and that one needs a non-perturbative fragmentation 
function (in addition to a charm distribution function in the photon). Since fragmentation 
and distribution functions are usually constructed in MS factorization, these terms must 
be subtracted also from the massive theory as soon as a fragmentation function, perturba- 
tive or non-perturbative, is taken into account. Therefore we studied the relation of this 
massive version with subtracted FSI terms with the massless 4-fiavour approach, using 
different assumptions concerning the fragmentation of the c quark into D* mesons. 

It turned out that the massive versions converged rather fast to their massless limits with 
increasing p-r- The convergence was strongest when using non-perturbative FF fitted to 
e~^e~ annihilation data 0. At low pt, for example at = 3 GeV, the massive cross 
section is reduced by approximately 20 % as compared to the massless approximation. At 
lower energies this reduction increases and makes the massless approximation unreliable. 
In addition the cross section is very much reduced by fragmentation effects, even at low 
Pt- This reduction increases from a factor of more than 3 at pt = 2 GeV to a factor of 
nearly 6 at = 12 GeV. Therefore, for reliable predictions one needs a good description 
of the fragmentation process. In our numerical evaluations we have taken the FF from fits 
to D* production in e~^e~ annihilation at LEPI. To compare with recent measurements of 
the inclusive D* cross section da/dpx in 77 reactions at LEPII we added the single- and 
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double-resolved cross sections, which up to now are available to us only in the massless 
approximation. The agreement of the calculations and the data is quite good down to 
Pt — 1-5 GeV even if we assume that the massless resolved contributions are reduced by 
a similar amount as the direct cross section. It remains to be seen whether a fit using 
additional experimental data to obtain a modified fragmentation function would lead to 
a better agreement. The fact that the data at medium pjn lie somewhat above the the- 
oretical predictions seem to indicate that a fragmentation function which is harder than 
the one we used in our analysis, could be favored. 

It is clear that in order to obtain results at even smaller pj- the finite charm mass cor- 
rections must be calculated also for the single- and double-resolved cross sections. Their 
knowledge would also be a prerequisite for the calculation of the total D* production cross 
section. 
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